Generalized Chaplygin Gas Model: 
Dark Energy - Dark Matter Unification and CMBR Constraints* 



M. C. Beutot, 0. Bertolami*, A. A. Sen§ 

Instituto Superior Tecnico, Departamento de Fisica, 
Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 



Abstract 

The generalized Chaplygin gas (GCG) model allows for an unified description of 
the recent accelerated expansion of the Universe and the evolution of energy density 
perturbations. This dark energy - dark matter unification is achieved through an 
exotic background fluid whose equation of state is given by p = —A/p a , where A is a 
positive constant and < a < 1. Stringent constraints on the model parameters can 
be obtained from recent WMAP and BOOMERanG bounds on the locations of the 
first few peaks and troughs of the Cosmic Microwave Background Radiation (CMBR) 
power spectrum as well as SNe la data. 



* Essay selected for an honorable mention by the Gravity Research Foundation, 2003. 

I 1 Also at CFIF, Instituto Superior Tecnico, Lisboa. Email address: bento@sirius.ist.utl.pt 

'■Also at CFNUL, Universidade de Lisboa. Email address: orfeu@cosmos.ist.utl.pt 

§Also at CENTRA, Instituto Superior Tecnico, Lisboa. Email address: anjan@x9.ist.utl.pt 



1 Introduction 



Cosmology is undergoing a blooming period. Precision measurements and highly 
predictive theories are coming together to yield a rich lore of data and methods 
that scrutinize existing models with increasing depth. It is quite remarkable that all 
available data can be fully harmonized within the Hot Big Bang Model, an unifying 
description in which several branches of physics meet to provide a consistent and 
testable scenario for the evolution of the Universe. In this picture, a particularly 
relevant role is played by Inflation, a period of accelerGiven the potential of the GCG 
model as a viable dark energy-dark matter unification scheme, many authors have 
studied constraints on the model parameters from observational data, particularly 
those arising from SNe la jTj and gravitational lensing statistics jH]. 

Quite stringent constraints arise also from the study of the position of the acoustic 
peaks and troughs of the CMBR power spectrum. The CMBR peaks arise from 
oscillations of the primeval plasma just before the Universe becomes transparent. 
Driving processes and the ensuing shifts on peak positions jHjated expansion in the 
very early Universe that allows for reconciling cosmology with causality and leads to 
a consistent explanation for the origin of the observed Large Scale Structure of the 
Universe. However, in order to fully account for the existing observations, one must 
bring in at least two additional new mysteries: the concept of Dark Matter, originally 
proposed to explain the rotation curves of galaxies and later used to address the issue 
of structure formation at large scales, and the idea of a smoothly distributed energy 
that cannot be identified with any form of matter, the so-called Dark Energy, needed 
to explain the recently observed accelerated expansion of the Universe. Even though 
these concepts are apparently unrelated, a scheme has emerged where an unification 
of these physical entities is possible through the rather exotic equation of state: 



where A a positive constant and a is a constant in the range < a < 1. This equation 
of state with a = 1 was first put foward in 1904 by the Russian physicist Chaplygin to 
describe adiabatic processes [T]; its generalization for a ^ 1 was originally proposed 
in Ref. [2] and the ensuing cosmology has been analysed in Ref. jHj- The idea that a 
cosmological model based on the Chaplygin gas could lead to the unification of dark 
energy and dark matter, thereby reducing two unknown physical entities into a single 



one was first advanced for the case a = 1 in Refs. jU E], and generalized to a 7^ 1 in 
Ref. 0. 

2 The Model 

The interesting behaviour of the equation of state (0) can be better appreciated 
by inserting it into the relativistic energy-momentum conservation equation, which 
implies for the evolution of the energy density [3] 

p*=( a+ j<m) ' (2) 

where a is the scale-factor of the Universe and B an integration constant. Remarkably, 
this model interpolates between a universe dominated by dust and a De Sitter one 
with an intermediate phase described by a mixture of vacuum energy density and a 
"soft" matter equation of state, p = ap (a ^ 1) 

Eq. (0) admits, in principle, a wider range of positive a values; however, the 
chosen range ensures that the sound velocity (c 2 s = aA/pl£ a ) does not exceed, in the 
"soft" equation of state phase, the velocity of light. Furthermore, as pointed out in 
Ref. [3], it is only for < a < 1 that the analysis of the evolution of energy density 
fluctuations is physically meaningful. 

More fundamentally, the model can be described, as discussed in Ref. jSj, by a 
complex scalar field whose action can be written as a generalized Born-Infeld action. 
This can be seen starting with the Lagrangian density for a massive complex scalar 
field, $, 

C = g^®„-Vm 2 ) , (3) 

which can be expressed in terms of its masss, m, as $ = (-^) exp(— im6). Assuming 
that the scale of the inhomogeneities is set by the spacetime variations of <fi corre- 
sponding to scales greater than m -1 , then <\>^ << m<f>, which, together with Eq.fT]). 
leads to a relationship between 2 and p: 

^{PcH)=p a ch {pT-A)^ , (4) 
and a Lagrangian density that has the form of a generalized Born-Infeld action: 



C-GBI = 



i - or 1+Q • (5) 



Notice that, for a = 1, one recovers the exact Born-Infeld action. It is easy to see 
that Eq. (J2J) has a bearing on the observed accelerated expansion of the Universe as it 
automatically leads to an asymptotic phase where the equation of state is dominated 
by a cosmological constant, 87r£L4 1//1+a , while at earlier times the energy density 
behaves as if dominated by non-relativistic matter. This dual behaviour is at the 
heart of the unification scheme provided by the GCG model. Figure 1 depicts the way 
the Universe evolves in the GCG model. It has also been shown that the underlying 
complex scalar field model admits, under conditions, an inhomogeneous generalization 
which can be regarded as a unification of dark matter and dark energy E] without 
conflict with standard structure formation scenarios [SI IU El El • It is clear that the 
GCG model collapses into the ACDM model when a = 0. 

These remarkable properties make the GCG model an interesting alternative to 
models where the accelerated expansion of the Universe arises from an uncancelled 
cosmological constant or a rolling scalar field as in quintessence models. 

In what follows, we shall discuss the observational bounds that can be set on the 
GCG model parameters. 



3 Observational Constraints 

Given the potential of the GCG model as a viable dark energy-dark matter unifica- 
tion scheme, many authors have studied constraints on the model parameters from 
observational data, particularly those arising from SNe la j7] and gravitational lensing 
statistics [H|. 

Quite stringent constraints arise also from the study of the position of the acoustic 
peaks and troughs of the CMBR power spectrum. The CMBR peaks arise from 
oscillations of the primeval plasma just before the Universe becomes transparent. 
Driving processes and the ensuing shifts on peak positions jU] 

Zpm =£A(rn- (p m ) , (6) 

where is the acoustic scale 



TO - Th 
= 7T — 

C s T h 



(7) 



t and Tig being the conformal time (r = / a~ l dt) today and at last scattering and c s 
the average sound speed before decoupling, are fairly independent of post recombi- 
nation physics and hence of the form of the potential and the nature of the late time 
acceleration mechanism. Hence, the rather accurate fitting formulae of Ref. ^0] can 
be used to compute the phase shifts <p m for the GCG model. In order to calculate 
the acoustic scale, we use Eq. (J2J) and write the Universe expansion rate as 



where A s = A/p]£o, Pcho = (A + B) 1 ' 1+a and we have included the contribution 
of radiation and baryons as these are not accounted for by the GCG equation of 
state. As discussed in Refs. [HJ[I2j, the above set of equations allow for obtaining 
the value of the fundamental acoustic scale by direct integration, using the fact that 



Comparing results from the above procedure with recent bounds on the location 
of the first two peaks and the first trough obtained by the WMAP collaboration |13j . 
namely £ Pl = 220 .1 ± 0.8, £ P2 = 546 ± 10, £ dl = 411.7 ± 3.5, together with the 
bound on the location of the third peak obtained by the BOOMERanG collaboration 
[HI, lp3 = 825^3, leads to quite strong constraints on the model parameters. These 
constraints can be summarized as follows |12j : 

1) The Chaplygin gas model, a = 1, is incompatible with the data and so are models 
with a > 0.6. 

2) For a = 0.6, consistency with data requires for the spectral tilt, n s > 0.97 and 
h < 0.68. 

3) The ACDM model barely fits the data for values of the spectral tilt 

n s ~ 1 (notice that WMAP data leads to n s = 0.99 ± 0.04) and for that h > 0.72 
is required. For low values of n s , ACDM is preferred to the GCG models whereas for 
intermediate values of n s , the GCG model is favoured only if a ~ 0.2. 

4) Our study of the peak locations in the (A s , a) plane shows that, varying h within 
the bounds h = 0.71±g;g| [ISj, does not lead to very relevant changes in the allowed 
regions, as compared to the value h=0.71 (see Fig. 3), even though these regions 
become slightly larger as they shift upwards for h < 0.71; the opposite trend is found 
for h > 0.71. 





= a 



-4 / da 




5) Our results are consistent with the bound found in Ref. jTTj using BOMERanG 
data for the third peak and Archeops data for the first peak as well as results 
from SNe la and age bounds, namely 0.81 < A s < 0.85 and 0.2 < a < 0.6. 

Bounds from SNe la data, which suggest that 0.6 < A s < 0.85 [ZJ, are also 
consistent with our results for n s = 1 and h = 0.71, which yield 0.78 < A s < 0.87. 

4 Discussion and Outlook 

In this essay, we have described the way the GCG model allows for a consistent de- 
scription of the accelerated expansion of the Universe and purports a scheme for the 
unification of dark energy and dark matter. This description is quite detailed and 
allows for an unambiguous confrontation with observational data. For this purpose, 
several studies were performed aiming to constrain the parameter space of the model 
using Supernovae data, the age of distant quasar sources, gravitational lensing statis- 
tics and the location of the first few peaks and troughs the CMBR power spectrum, 
as measured by the WMAP and BOOMERanG collaborations. These studies reveal 
that a sizeable portion of the parameter space of the GCG model is excluded. 

More concretely, our results indicate that the Chaplygin gas model, a = 1, is 
incompatible with the data and so are models with a > 0.6. For a = 0.6, consistency 
with observations requires that n s > 0.97. We find that the ACDM model hardly 
fits the data for n s ~ 1 and h > 0.72 is required. For lower values of n s , ACDM is 
preferred to the GCG models whereas for intermediate values of n s the GCG model 
is favoured only if a ~ 0.2. 

We conclude that the GCG is a viable dark matter - dark energy model in that 
it is compatible with standard structure formation scenarios. Moreover, although its 
parameter space is rather constrained, the model is consistent with all the available 
Supernovae, gravitational lensing and CMBR data. Finally, the model does not suffer 
from the well-known fine-tuning problems that are present in alternative dark energy 
candidate theories such as ACDM and quintessence models. 
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Figure 1: Cosmological evolution of the Universe described by the Generalized Chap- 
lygin Gas model. 




Figure 2: Contour plots of the locations of the first three peaks and first trough of 
the CMBR power spectrum, in the (A s ,a) plane, for a GCG model, with h = 0.71, 
for different values of n s . Full, dashed, dot-dashed and dotted contours correspond 
to observational bounds on £ pi , £ P2 , £ P3 and i^, respectively. 



